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Abstract: We initiate a systematic study of entanglement and Renyi entropies in the 
presence of gauge and gravitational anomalies in even-dimensional quantum field theories. 
We argue that the mixed and gravitational anomalies are sensitive to boosts and obtain a 
closed form expression for their behavior under such transformations. Explicit constructions 
exhibiting the dependence of entanglement entropy on boosts is provided for theories on 
spacetimes with non-trivial magnetic fluxes and (or) non-vanishing Pontryagin classes. 
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1 Introduction and summary 

A remarkable feature of quantum field theory is that an admissible classical symmetry may 
be broken by quantum mechanical effects. When a theory possesses such a broken symmetry 
we refer to that symmetry as anomalous. 

There is a substantial volume of scientific works on the effect of anomalies on S-matrix 
elements, see, e.g., [1-3], for several reviews. However, little is known regarding the mani¬ 
festation of anomalies in other field theoretic quantities. Recently, it has been realised that 
anomalies play a significant role in the thermodynamic behaviour of anomalous gauge theories 
[4-11] leading to novel effects in heavy ion collisions (see, e.g., [12-14]), specialized Weyl semi¬ 
metals [15-18] and astrophysical phenomenon [19-21]. In this work, we focus our attention 
on the effect of anomalies on entanglement entropy, previously studied also in [22-25] . 

If the Hilbert space "H of a quantum mechanical system can be decomposed into a direct 
product % = then the entanglement entropy Sa of a state \'ijj) provides a measure of 

correlation between degrees of freedom of 1^) in T-La and in T-La- More formally, to compute 
the entanglement entropy one traces over the degrees of freedom in the complement A to 
obtain a density matrix pA = tr^|V’)(V’| in the subsystem A. The entanglement entropy for 
the region A is given by the von Neumann entropy associated with pA, 

S'a =-tr^ (pa log pa) ■ (1-1) 
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In quantum field theory one often considers the entanglement entropy associated with 
a spatial region A. Starting with a state \'ip) the entanglement entropy is given by the von 
Neumann entropy of the reduced density matrix pA obtained by tracing over the Hilbert space 
of the spatial region A. In what follows we will refer to the interface between A and A as 
the entangling region S = dA = dA and work under the assumption that the Hilbert space 
of the theory can be decomposed into a tensor product of the Hilbert space associated with 
A and that associated with A. We will refer to (the Euclidean version of) the manifold on 
which the theory is placed as M. 

In order to compute the entanglement entropy Sa in a quantum held theory one uses 
(1.1) on a regulated version of the theory. When the initial state \'ip) is the vacuum state |0), 
the continuum limit of (1.1) is given by 


Sa = lim Sn , 

n—>■! 


(1.2a) 


where the Renyi entropy, Sn, is given by (the analytic continuation of) 


Sn = - 


nWi - Wn 
n — 1 


(1.2b) 


For integer n, Wn = — In D(j)e~^ is the Euclidean generating function on the manifold 
Ain where Ain is the n-fold cover of Ai such that each sheet of the cover is connected 
along A [26]. Therefore, Ain has a codimension-two singularity at the entangling surface T, 
with surplus angle 27r(re — 1). For non-integer n, Wn is dehned by analytically continuing 
the integer valued Wn- Note that Ad = is the original Euclidean space without any 
singularity. The analytic continuation from integer n > 2 to n = 1 used in (1.2a) is usually 
carried out under an implicit assumption regarding a replica symmetry, Zn, between sheets 
in Ain- See e.g., the discussions in [27-29] for the role of the symmetry in the holographic 
formula of entanglement entropy [30, 31]. In the remainder of this work we will assume that 
the analytic continuation from non-integer n to integer n can be carried out and take (1.2a) 
to be the definition of the entanglement entropy in quantum field theory. Our conventions 
for Wick rotation are summarized in appendix A. 

From a formal perspective the definition (1.2a) implies that entanglement entropy is 
associated with the Euclidean partition function of the theory in the limit where the manifold 
on which the theory is dehned becomes regular. It is then perhaps not surprising that the 
entanglement entropy is sensitive to anomalies. 

Indeed, let us consider a theory with gauge, gravitational, and mixed gauge-gravitational 
anomalies. By gauge and gravitational anomalies we mean anomalies associated with turning 
on a background (non-dynamical) source term for the current and stress-energy tensor. In 
what follows we will refer to these source terms as a gauge-held and metric respectively, but 
we emphasize that these are non-dynamical helds and the associated anomalous symmetries 
are global. Anomalies associated with local symmetries generate non-unitary theories at best 
and will not be discussed further. For simplicity, we focus on Abelian gauge anomalies but 
our results can be easily generalized to non-Abelian ones. 
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In the presence of gravitational or mixed anomalies the stress-energy tensor is no 
longer conserved. Instead, one finds the non-conservation law 




(1.3) 


where is a field strength associated with possible external fields, is the current as¬ 
sociated with a global Abelian symmetry and denotes the anomalous contribution to the 
(non-)conservation law of the stress-energy tensor. The term is associated with Joule 

heating. The explicit form of can be computed using the Wess-Zumino consistency condi¬ 
tions. It is given by 

{^gr%) , (1.4) 

where is the gauge field associated with F^y and J and may be expressed as variations 
of the Chern-Simons forms associated with the anomaly. In the presence of boundaries (1.4) 
is expected to receive corrections [32-34]. We provide a derivation of (1.4) in section 2. A 
derivation of the explicit expressions for J and can be found in appendix B and in [10]. 
It is worth emphasizing that in (1.4) specifies the non-conservation law for the consistent 
stress tensor in Euclidean signature. The expression for r associated with the covariant stress 
tensor is somewhat different. 

Consider boosting the state |0) by a hyperbolic angle k which implies a rotation of the 
Euclidean vacuum by an angle 9 = in. Let be the generator of this rotation. We argue 
that if the dynamics are invariant under 9 in the absence of an anomaly, then the Renyi 
entropy associated with the boosted state |0) in the presence of the anomaly satisfies 


deSn 


e=o 


. J Ml J Mn 


(1.5) 


Equation (1.5) is valid in the absence of boundaries. We discuss possible subtle corrections 
to (1.5) for manifolds with boundaries or coordinate singularities in section 3. A similar 
expression for the entanglement entropy can be obtained by inserting (1.5) into (1.2a). In 
what follows we will omit the symbol, keeping in mind that all our results for derivatives 
of the entanglement entropy are evaluated at d = 0. The imaginary contribution to the 
entanglement entropy follows from the analytic continuation k —)• —i9 and is removed when 
analytically continuing back to the Lorentzian signature manifold. 

We carry out an explicit evaluation of (1.5) in configurations where the entangling surface 
splits space into two regions, i.e., A = {x\x^ > 0}. In this case the (Euclidean) metric on A4n 
is given by ds^ = dp^ J- p^dr'^ + gij{x) dx'^dx^ where 0 < r < 27rn. Such a geometry 

describes a cone with opening angle Inn. Using (1.5) we obtain the following results: for a 
two-dimensional quantum field theory and the vacuum state we have 


leading to 


deSn = 2-Kicg 



deSA = 47rf Cg , 


( 1 . 6 ) 

(1.7) 
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where Cg = — Xli Xi/(967r) is the strength of the anomaly such that the index i runs over 
all chiral particle species with Xi denoting the chirality. For a conformal field theory Cg = 
[cl — cr) /(967r). In order to make contact with the existing literature we argue, in section 2.1, 
that if we use a finite interval instead of a semi-infinite line for the entangling surface then 
(1.7) will be modified to dgSA = Snicg. A comparison between our field theory results and 
a holographic computation can be found in appendix B.3. A more elaborate discussion and 
comparison can be found in section 3 where we discuss possible corrections due to boundary 
terms. 

For a four-dimensional theory and in the presence of an external magnetic field, F 23 = B, 
through the entangling surface, we obtain 

dgSA = dvri a CmB Vols , (1-8) 

where Cm = Xi9i/(1927r^) is the strength of the anomaly (with qi denoting the fermion 

charge and Xi the chirality) and ol denotes an ambiguity in the Chern-Simons term derived 
from the mixed anomaly polynomial F f\ R f\ R and Vols is the volume of the entangling 
surface S. If we choose the contact terms in the theory such that the anomaly is manifest 
only in the gravitational sector, then a = 1. 

For a six-dimensional theory in Minkowski space and with a magnetic field F 23 = Bi and 
F 45 = B 2 , we obtain 

dgSA = dTTZ a CmBiB 2 Vols , (1.9) 

where Cm = Xj 9 f/( 7687 r^) is the strength of the mixed anomaly in six dimensions (with 

Xi and Qi denoting the chirality and charge of fermion species). If we place the theory on 
X i^3, we find 

dgSA = Ca , ( 1 - 10 ) 

where now Ca = — YliiXi ~ 8ti)/(368647r^) (with Xi and ti counting fermion species chirality 
and self-dual and anti-self-dual two-form field species) is the strength of one of the two 
possible gravitational anomalies of a six-dimensional theory. It seems that in order to obtain 
a non-zero value for the entanglement entropy of the other gravitational anomaly one needs 
to consider more intricate geometries. In section 2.4 we present results for higher-dimensional 
theories. 

While this work was being completed we became aware of [35] where it was pointed out 
that coordinate singularities may become physical once gravitational anomalies are present. 
We thank the authors of [35] for sharing their draft with us prior to publication. We discuss 
the possible role of such corrections in section 3. 

2 Entanglement and anomalies 

Equation (1.2a) implies that the entanglement entropy is associated with the limiting be¬ 
haviour of the generating function as the manifold it is placed on becomes non-singular. 
Indeed, the behaviour of entanglement entropy on scaling transformations of the entangling 
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surface can be tied to the scaling behaviour of the generating function Wn which, in turn, de¬ 
pends on the central charges of the theory if the latter is conformal [31, 36-39]. In particular, 
if we denote a scaling variation of the entangling surface by 5^ then 


5aSA = — lim 

n—^1 


n5^Wi - 6^Wn 
n — 1 


(2.1a) 


The scaling is related to the trace of the stress tensor, 

SaWn = \ j , 


(2.1b) 


which is related to the central charges, Cj and a, of the theory when the latter is conformally 
invariant, 

T% = 2i-l)'^/^aE-Y,c^Ii. (2.1c) 

i 

By evaluating the curvature terms Ij and £" on a spherical entangling surface one can infer 
the dependence of the entanglement entropy on the central charges directly from (2.1). See 
[31, 37, 38, 40] for details. 

Likewise, we can consider the behaviour of the generating function under a coordinate 
transformation which generates an isometry, e.g., a boost or rotation. If we denote the 
generator of boosts by 6o then 


deSn = 


n6eWi - SgWn 
n — 1 


( 2 . 2 ) 


The behavior of Wn under boosts (and coordinate transformations in general) is fixed by 
the Wess-Zumino consistency condition. In appendix B we argue that for a small coordinate 
transformation ^ one has 


= i 


d<^x^d,x^T%, 


(2.3) 


with 




dies 


(2.4) 


where Ics is the (Euclidean) Chern-Simons form associated with the anomaly, is related 
to the Christoffel connection via T^i, = T^undx^ and ^'T^u is the volume form dual to 
The derivative d/dV is taken at constant R and picks up a minus sign when acting on 
p-forms similar to the exterior derivative d. 

The energy-momentum tensor and conserved current can be obtained by varying the 
generating function with respect to the metric and gauge field respectively. Thus, under a 
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small coordinate transformation parameterized by ''^6 have 


< 5 ;, 144 . = j > 

= I [V^Xu {X'^F,^ + d^ix'^A,))] , (2.5) 

= - I d'^x^x. (V^T^" - , 


where in the last step we have integrated by parts. Comparing (2.5) to (2.3) yields (1.4) 
where we have defined 


V;. = -ij . 


( 2 . 6 ) 


A non-zero J implies that the theory is not invariant under gauge transformations. Indeed, 
carrying out an analysis similar to the one carried out in appendix B but for gauge transfor¬ 
mations results in 


*J 


dies 

dA^ ■ 


(2.7) 


See, e.g., [10] for a full derivation. 

We note in passing that in the presence of boundaries, (2.3) and (2.5) receive corrections 
leading to a modification of our main result. On one hand, boundary terms generated by 
integrating (2.5) by parts will contribute to the conservation law in addition to a possible 
conservation law for the boundary stress tensor. On the other hand, a boundary term may 
appear in the right hand side of (2.3) due to the anomaly. Such a boundary term may be 
determined by solving the Wess-Zumino consistency condition. 

In constructing the generating function 144i one must deal with ambiguities associated 
with possible contact terms (often referred to as Bardeen counterterms [41]) which can shift 
the appearance of the anomaly in non-conservation laws. For instance, in the presence of a 
mixed gauge-gravitational anomaly, one can choose an appropriate Bardeen counterterm so 
that the anomaly will manifest itself entirely in the gravitational sector. A more familiar setup 
where a Bardeen counterterm is present is the (mixed) axial-vector anomaly; the existence of 
the Bardeen counterterm is what allows us to enforce the current conservation in the vector 
sector. In the current context the ambiguities relevant to our computation are captured by 
a choice of the Chern-Simons form Iqs- We will provide explicit examples of such in the 
following subsections. An extensive discussion of these terms in the language of the current 
work can be found in, for example, [42]. 

We now choose x^ to specify a boost by a (small) constant hyperbolic angle 56 in the r 
rapidity direction. We would like to compare the entanglement and Renyi entropies in boosted 
and unboosted entangling regions. That is, we would like to compute, say, the Renyi entropy 
for an unboosted entangling region and compare it to the Renyi entropy for a boosted one 
while keeping the cutoff parameters fixed. When constructing the n-fold cover of the boosted 
entangling surface, each sheet of the covering will be shifted by an angle 56 such that, overall. 
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Airi rotates by an angle n66 (implying that rotating by a Euclidean angle 9 = 2'k will bring 
us back to the original state). Thus, 

dW 

= W{e + n59) - W{e) = nse— . (2.8) 

o9 

In what follows, we will denote = n59£^^. Inserting (2.8) into (2.3) and using (2.2) we 
obtain (1.5). 

The careful reader will note that in the coordinate system in which = dr one 

has = 0. For this reason we will work exclusively in a Cartesian coordinate system 

where d^^^ ^ 0. The lack of invariance under coordinate transformations is a hallmark of 
gravitational anomalies. Indeed, the non-tensorial behavior of implies that the stress 
tensor derived from the generating function Wn (referred to as the consistent stress tensor) 
does not transform like a tensor under coordinate transformations. It is possible to define a 
covariant stress tensor by adding to the consistent stress tensor a polynomial which is local 
in the external fields. Such a stress tensor however can not be derived from a local generating 
function and therefore its role in entanglement entropy is somewhat obscure. 

In the remainder of this work, we evaluate the change in the R&yi entropies, (1.5), for 
geometries of the form x AA, and the particular case of an entangling region which divides 
space into two equal halves, A = {x\x^ > 0} where the x^ direction is along the spatial part 
of the Lorentzian continuation of into 

ds^ = dt^ + (dx^)^ + dsj^f , (2.9) 

where dsfj- = metric of the {d — 2)-dimensional manifold parame¬ 

terized by the coordinates Xi (i = 2, 3 • • • , d — 1). We will consider slightly more complicated 
geometries in a future publication. The metric (2.9) describes a cone of opening angle 27rn 
and can be rewritten in the form 


ds — dp p dx dsjij ', 


( 2 . 10 ) 


where r has periodicity 27rn. The (Euclidean) cone located at p = 0 needs to be regulated 
in order to carry out explicit computations. In what follows we will smoothen the tip of the 
cone by using a regulating function u such that the metric takes the form [43] 

ds^ = u{p)dp^ + p^dr"^ + ds\ , (2-11) 


where tt(0) = and u asymptotes to unity at large p outside a small region of radius e 
such that u{p) = 1 and u'{p) = 0 for p > e. See [43] for details. Going back to Cartesian 
coordinates we find that the regularized cone metric takes the form 

( 2 . 12 ) 


,n [tdx^ — x^dtY + u{t,x^){tdt + x^dx^y ,n 
ds^ = - '—7. —V-tttP- — + dsjj-. 


+ {x^y 

Note that the origin t = 0 is not covered by either (2.11) or (2.12). Had we used the 
coordinate system advocated by [27] 


ds^ = -|- (x^)^ + a^y (dt^ + d(x^)^) , 


(2.13) 
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with e = re — 1 and a regularization parameter a, as we do in appendix B.2 then there is no 
coordinate singularity at the origin. 

Following our conventions for Wick rotation, (A.9), the generator of rotations 
(boosts analytically continued to Euclidean signature) can be related to via 

X^5^ = -n 66 {tdx - xdt) = , (2.14) 


implying 



(2.15) 


Inserting (2.15) into (2.3) we find that 

5eW^ = inSe j dtdx^ J d^-^Xy/^ (-T^xi + • (2-16) 

A further simplification of (1.5) can be obtained by noting that on A4i there is no conical 
singularity and therefore T^v vanishes there, 


[ d^x^d,eT\ = Q. 

JMi 


(2.17) 


Inserting (2.16) and (2.17) into (1.5) we find that, in the geometries we are considering. 


deSnl-o = (- 7^.1 + r^\) • 

n-lJMr, ^ > 


(2.18) 


We now proceed to compute (2.18) in various dimensions. 


2.1 Two-dimensional theories 


For two-dimensional theories there is a single type of gravitational anomaly whose Chern- 
Simons term I 3 can be derived from the four-dimensional anomaly polynomial 


Cg R^u A R!' g = Cgd 


F'^^ A R% - A V% A TPg 

O 


= dh, 


(2.19) 


where RP^ = ^RP^padxP A dx^ is the Riemann tensor two-form and F^jy = T^ypdxP is 
the Christoffel one-form. The two-forms are related by R^y = dT^y + T^p A T^y. The 
multiplicative constant Cg can be determined explicitly from the field content of the theory, 


where the sum is over all species and Xi is the chirality of the fermions and scalars of the 
two-dimensional theory (Majorana-Weyl fermions contribute ±1/2). For a conformal field 
theory one has 


CL - cr 






where cl and cr are the left and right central charges respectively. In what follows we find it 
convenient to omit the wedge product and replace the spacetime indices with a trace so that 
(2.19) takes the form 


Tr (i?2) = cg dTr ( ri? - ) = dh • 


The explicit form of T^u can be obtained from (2.4) as 


d 


1 , 


= Tr(^ri?--r^j =c,dr, 

resulting in 

Evaluating (2.15) and (2.24) on the regularized cone (2.11) one finds 

J y/gT^vdtdx^ = — J ^^^^dtdx^ , 

where = t^ + (x^)^ and we have used y/g. Using (2.16) this yields 


deWn = 27r(n^ — l)id0Cg 


( 2 . 22 ) 

(2.23) 

(2.24) 

(2.25) 

(2.26) 


The complex value of 5gWn is expected given that the analytic continuation of —i69 back to 
Lorentzian signature is given by 

deWn —)• —27r(n^ — l)SKCg . (2.27) 

where 5k is an infinitesimal boost parameter. 

Inserting (2.25) into (2.18), we obtain 

deSn = 2Tricg(^l + ^^ , (2.28) 

which implies 

dgSA = dvri Cg , (2.29) 

as argued in (1.7). 

We note that it is also possible to use our formalism to compute the change in the 
entanglement entropy on a finite interval of length L. To this end, consider a manifold Ain 
with two conical singularities such that the metric takes the form 

ds'^ = dt^ + (dx^)^ , (2.30) 


away from + (x^)^ < e and away from + (x^ — T)^ < e, 

(tdx^ — x^dtY + u + (x^)^^ (tdt + x^dx^)'^ 


ds^ = 


+ {x^Y 


(2.31) 
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for ^ 


ds^ = 


{tdx^ — — L)dtY + u + {x^ — L)'^^ (tdt + (x^ — L)dx^Y 

+ (x^ — Ly 


(2.32) 


for y + (x^ — Vy < £. The line elements for (2.31) and (2.32) have been obtained from (2.11) 
with appropriate shifts. Since T^u vanishes for y + (x^)^ > s and ^ + (x^ — L)^ > e we need 
to evaluate (2.16) only near the conical singularities. A quick computation yields 


dgSn = dvricg ( 1 + 


1 


n 


which implies 


(2.33) 

(2.34) 


QqSa = Svri Cg . 

In appendix B.3 we compare the results of this section to a holographic computation. 

2.2 Four-dimensional theories 

Four-dimensional theories do not have a gravitational anomaly but may have a mixed gauge- 
gravitational anomaly. Thus, we may write 


,FTr(i?2) =cmd aFTi{TR-^T^) + {1 - a) A Tv {R^) 


= dh 


(2.35) 


The parameter a is a free parameter which manifests the possibility of adding a (Bardeen) 
contact term to the generating function and Cm is the strength of the mixed anomaly given 
by 

( 2 - 36 ) 


Cm. — 


with Xi 3.nd qi the chirality and charge respectively of the fermion species in the theory. As 
we will see shortly, such a contact term offers the possibility of shifting the anomaly entirely 
into the gauge sector of the theory. Indeed, using (1.4), one has 




implying a (non-)conservation law for the (Euclidean) stress tensor of the form 


(2.37) 


+ iJA^^ + iVgTPy • 

yg 


(2.38) 


The U{1) current coupled to the external source satisfies the (non-)conservation law 

= (2.39) 
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When a = 1 the current is conserved and the anomaly is manifest only in the energy- 
momentum tensor. 

If we are to use (1-5) then we need that T^u 7 ^ 0. To this end we consider the entanglement 
entropy associated with a state |0 )b where B = T 23 is a magnetic field in the direction. 
In this case, 

J dx'^dx^ = a , (2.40) 

with <I> the magnetic flux through M at x^ = 0, i.e., the entangling surface. Denoting the 
volume of the entangling surface by Vols, the flux can be written in the form $ = i?Vols. 
Expression (2.40) is identical to (2.24) upon identifying Cg = acm^- Thus, 

deSn = 2'KiacmB Vols , (2-41) 

leading to 

SoSa = 47 rz a CmB Vols • (2.42) 

The factor of a in (2.41) implies that the entanglement entropy depends on the particular 
Bardeen counterterm used to shift the anomaly between the gravitational sector and the gauge 
sector. Indeed, we expect that if the gravitational anomaly does not appear in the appropriate 
part of the non-conservation law for the stress tensor then the entanglement entropy would 
not be sensitive to it. As discussed earlier, the particular choice of a in a given theory depends 
on how the currents are coupled to dynamical gauge fields or to the metric. Much like the 
canonical axial-vector anomaly, if we plan on coupling to a dynamical gauge field then we 
will be forced to set a = 1 . 


2.3 Six-dimensional theories 


Six-dimensional theories have two types of gravitational anomalies and a mixed gauge-gravitational 
anomaly. We classify the Chern-Simons terms associated with these anomalies by 


Tr (i^2) =cmd (1 - a)AF Tr {R^) + oF^Tr ( TR - 

O 


= dF 


7,m 


Ti{R^)^ = Cad 


Tr ( TF - -r^ ) Tr {R^) 


Cb Tr (R‘^) = CbdTr 


R^r - -R^r^ - -RF^Rr + -rf^ - 


— dll^a ) 
1 


35 


= dL 


7,b ■ 


Here 


Cm. — 


TGStt^ 


Xiqf , 


1 


Ca = - 


SGSGdvr^ 


Y ’ 


(2.43) 


Cb = 


^J2iXi + 28ti) , 
(2.44) 


4G0807r3 


where Xi s-nd qt denote the chirality and charge of the fermion species and ti counts the 
number of self-dual (and anti-self-dual) two-form fields. 
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Using (1.4) we find that 


= CmiTmTu + CaiTaYu + C^TbYv , 


(2.45a) 


where 

{TmYu = F.sF^r,, 

{TaYu = 


pnr) 1 


{%Yu = I 


filKri 


7 

_DM 2 cT/^S pM 4 

-„4X K/X 4 J- TJlll 


(2.45b) 


10 


I r/^2 pM3 J pM4 pMS 
W J- 7/43^- (5/^4 KMS^ V^^l 


Let us consider a manifold (2.9) which is of the form x Af. Since the Riemann tensor 
two-form and the Christoffel connection one-form are going to be block diagonal then {TbYv 
given in (2.45b) will vanish. Taking advantage of the lesson learned from the analysis of the 
mixed anomaly in four dimensions we consider a manifold (2.9) with a constant magnetic flux 
through M, i.e., 

/ = $ . (2.46) 

Jm 

Then, integrating F^u over M we obtain an expression similar to the two-dimensional result 
(2.24) but with Cg replaced by acm^ — 247r^Car[S], where r[S] is the Hirzebruch signature on 
the entangling surface, 

Thus, 


deSn = 27ri 1 -|- 


^1 -F (acm$ - 247r^Car[E]) , 


(2.48) 


and 


deSA = Fki {acm^ - 247r^Car[S]) . (2.49) 

To clarify the role of and r[S], let us consider two possible setups. If AA = T^, 
F 23 = —F ^2 = Bii A 45 = —F 54 = B 2 and the remaining F's vanishing, we obtain r[S] = 0 
and <I> = Riil 2 Vols with Vols the volume of the entangling surface. Likewise, if F = 0 and 
J\f = F3 then r[S] = —16 and $ = 0. 

It seems that in order to observe the effect of the second type of gravitational anomaly the 
topology of the n-fold cover of M. must not reduce to a (direct) product manifold. While it 
is straightforward to conjure a manifold for which the last line of (2.45b) does not vanish, its 
interpretation as an n-fold cover of a manifold A4 and therefore its relation to entanglement 
entropy is somewhat obscure. We discuss this case further in section 3. 
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2.4 Higher dimensions 


We are now in a position to compute the entanglement entropy associated with an even 
d-dimensional theory. The anomaly polynomial for the theories we are considering will be 
characterized by a set of integers ki and and take the form of 

^ c, J] Tr = did+i , (2.50) 

i n=l 


where for each i, 2ki + = d -|- 2. The analysis of the previous sections implies that 

for a manifold of the form x Af, we will have 


where 


and 


deSn = 27ri ( 1 -|— 
\ n 

)a. 

(2.51) 

C = ^ Ci [ F^^ Tr(i72)™i- 

^ Yl Tr(i72^)"*- , 

(2.52) 

n=2 


1 Cl 5 

CD 

II 

(2.53) 

Ci = { 


[UiCi (no sum), 

ki> 0, 



where the at specifies a free parameter which determines in which sector the mixed anomaly 
will be manifest. Equation (2.51) leads to 


deSA = dvriC. 


(2.54) 


One can single out the contribution of a particular coefficient Cj with > 0 in (2.50) by 
an appropriate choice of the manifold Af. Since the presence of anomalies is inherently related 
to chiral fermions (or self-dual p-form helds) it would be convenient to choose an Af which 
can be endowed with a spin (or string) structure.^ For instance, choosing a spin manifold 


Af = X X (HP’^)”"" , (2.55) 

n=2 

implies that Sn and Sa will be proportional to cj only. Using our former result for a K3 
surface 



Tr (i?2) = 96(27r)2 , 


and 


[ Tr(i?2^) =2(2n-b2-4”)(-47r2)”, 

JHIP" 


for quaternion projective spaces [44], we find that 


C = cj^j (96(27r)2)'"i ^ Yl (2(2n + 2 - 4^)(-47r2)^)' 


n=2 


(2.56) 


(2.57) 


(2.58) 


with 

^We thank K. Jensen and Y. Tachikawa for discussions on this point. 
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3 Discussion 


The entanglement entropy associated with a state |?/>) and entangling region S in a quantum 
field theory is dehned by (1.1) where pA is the reduced density matrix computed at a particular 
instant of time, t = to- If the state \'ip) is an eigenvalue of the Hamiltonian H then its time 
evolution is trivial then the entanglement entropy will be independent of to- Likewise, if \'ip) 
transforms trivially under boosts, one expects that Sa will be invariant under boosts as well. 
More precisely, we expect that two inertial observers using identical entangling surfaces and 
cutoff schemes will agree on the entanglement entropy Sa- 

To make a sharper statement about the boost invariance of the entanglement entropy let 
us take a closer look at the comparison of the entanglement entropies between two inertial 
observers, call them a and b. Let observer b move at constant rapidity k relative to observer 
a. Both observers are interested in computing the entanglement entropy associated with an 
entangling surface S and both observers use the same cutoff scheme. Suppose observer a 
computes the entanglement entropy 5 'a( 0) in his own frame and the expected entanglement 
entropy S'a(k) in the frame of b. 

In his own frame, S'a(O) can be computed using the replica trick (1.2a), that is, by 
computing the partition function on Mn, the n-fold cover of Ai where each of the sheets 
of Ain are connected along the entangling region. According to observer a, Sa{h) can be 
computed by considering the partition function on Ain(fi), the n-fold cover of Ai where each 
of the sheets of Ain{G) are connected along an entangling region which is rotated by the angle 
6 = —in relative to the original. 

As we have seen, in the presence of anomalies one hnds that S'a(^) and ^^(O) are not 
necessarily equal. In particular, we have seen that for (Euclidean) manifolds of the form 

= dt^ + (dx^)^ + ds^f , (3.1) 

and an entangling surface S = {x 11 = 0, = 0}) th® change in entanglement entropy due to 

a boost is given by 

dgSA = dvri (7 , (3.2) 

where C is given by (2.58) and depends on the details of the anomaly and the manifold M. 
We have argued that one may choose a manifold N such that C will be non-zero as long as 
the anomaly polynomial P satisfies, 

i.e., it’s derivative with respect to Tr with all other Tr held fixed is non-zero. 

It is unclear whether the entanglement entropy is susceptible to anomalies for which the 
right hand side of (3.3) vanishes. From the arguments presented in section 2.3 we expect that 
for the latter type of anomalies one would need to consider more intricate entangling surfaces 
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S, or manifolds A4 whose structure is different from (3.1). For instance, if we use 


ds^ = 


dp^ + {dr + uji{y^)dy^ + uj 2 {y^)dy'^Y + '^{dy" 


(3.4) 


2=1 


as the line element for a six-dimensional manifold Ain with 0 < r < 27rn then it is a 
straightforward (though tedious) exercise to show that at least for small a;*, dgSA will receive 
contributions from both Ca and Ch defined in (2.45), the latter being associated with a Tr(i?^) 
term in the anomaly polynomial. It would be interesting to be able to associate a well-defined 
entangling surface with (3.4) or with a variant of it. 

A result similar to that in (1.5) appeared in [23] for two-dimensional conformal field theo¬ 
ries with a finite entangling surface. By extending the standard argument relating correlation 
functions of twist operators to the partition function on Ain, Zn = [26, 45], it was 

argued in [23] that 


Zn — 


Cn 


^2h^2h ’ 


with 


^ CL f 1 

24 V n 


z CR f 1 

24 V n 


(3.5) 


(3.6) 


Here the entangling surface consists of the points 0 and z in Ain and Cn is a constant which 
does not depend on z. Let us consider z = R and carry out an infinitesimal coordinate 
transformation z R{1 + ind9). On one hand, following the general analysis of section 2.1,^ 
the change in the partition function due to rotations by an angle n56 is given by 


SgZn 


= —i56{n^ — 1 ) 


CL - Cr 

24 


(3.7) 


On the other hand, we have 


^2h^2h ) ^ ^ SCn 


Cn 


= —i59{n‘^ — 1 )- 


12 


+ 


Cn 


(3.8) 


The validity of (3.5) seems to imply that 

SCn 


Cn 


= i59{n^ — 1 ) 


Cl - Cr 

24 


(3.9) 


A computation of correlation functions for twist operators in the absence of a gravitational 
anomaly has been carried out in [46] where an explicit expression for Cn has been obtained. 


^Equation (3.7) can be obtained by using the same manipulations as those which lead to (2.34). Note that 
the coordinate transformation {t, x^) —^ {z, z) with z = +it and z = x^ — it is a linear orientation reversing 
coordinate transformation so that = ~^l\f§ relative to = 1/y^. Here is the square root of the 
determinant of the metric which is imaginary in the z, z coordinate system and real in the t, x^ coordinate 
system. 
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By adding the contribution of the effective action for gravitational anomalies [41, 47, 48] to 
the analysis of [46] one should be able to verify the form (3.5) and derive (3.9). 

Note that if we set 6 Cn = 0 then (3.8) will differ from (3.7) by a factor of 2. Curiously, 
such a factor of 2 will appear in the holographic computation of the entanglement entropy if 
we mistreat certain terms which are integrated by parts (see the discussion in appendix B.2). 
Further, we can define a covariant entanglement entropy by using 

6 gWi^°^^=i[ (3-10) 

JMn 

where 

= Cg e^'^dyR , (3-11) 

and R is the Ricci scalar. The expression (3.11) follows from the non-conservation law for 
the covariant stress tensor, 

V^r(c°v)/^'^ =(3.12) 
in the absence of an external field strength F. We find that 

= 47r(n2 - 1) m Cg , (3.13) 

which is twice as large as (2.26). Inserting (3.13) and 6 gWi = 0 into (1.2b) we obtain 

= 47rz (^1 + ^Cg. (3.14) 

A possible resolution of the aforementioned factor of two discrepancy between the results 
presented here (and the holographic result) and the conformal field theory computation pre¬ 
sented above has been advocated in [35]. Since the origin of our coordinate system t = = 0 

contains a coordinate singularity it is possible that the origin contains an extra delta function 
contribution. While such a contribution is, perhaps, surprising given that the geometry has 
been regulated we can not rule out such a possibility. As we pointed out earlier, one can 
construct a coordinate system such as (2.13) which can be used to compute the entanglement 
entropy and does not posses a coordinate singularity. Such a coordinate system may a priori 
resolve the problems raised by [35]. However, in [49] it was argued that one would expect, 
on physical grounds, that the entangling region be smeared over an ultraviolet parameter 
£ which, in essence, serves as a boundary. It may be that such arguments, when applied 
to anomalies, introduce contributions to the entanglement entropy which persist even as we 
take the limit where the boundary vanishes. A preliminary analysis of the modification of 
entanglement entropy due to boosts in the presence of boundaries has been carried out in [35] 
where some, but not all, contributions of boundary terms have been analyzed. We postpone 
a full analysis of boundary terms for future work. 

Our result (1.5) implies that the entanglement entropy is sensitive to gravitational and 
mixed gauge-gravitational anomalies. It would be interesting to understand whether the en¬ 
tanglement entropy is also sensitive to discrete anomalies such as parity, or to gauge anomalies 
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(or gauge transformations in the presence of mixed anomalies). Indeed, following an anal¬ 
ysis similar to the one we have used so far, one expects that under a gauge transformation 
A ^ A + d\ the generating function transform as 


6xWn = 


\*J. 


'Mr. 


(3.15) 


For a 17(1) polygon anomaly in d 
to 


2m spacetime dimensions and A a constant, (3.15) reduces 


5xWn = icx[ F^\ 
JMn 


(3.16) 


where c is the strength of the anomaly. 

Since we have defined F to be the (external) field strength on A^i, it is insensitive to 
the n-sheeted cover and then 5xW oc n in which case dxSn = 0. Put differently, the external 
gauge field A has the replica symmetry A(r = 0) = A{t = 27r) = • • • = A(r = 27rn) 
on Ain which yields F^ = n F"^ for an integer n. Inserting the former result into 
(3.16) and using (1.2b) results in dxSn = 0 for integer n > 2 implying OxSa = 0. 

A similar argument, relying also on the behavior of curvature invariants on Ain [43] 
implies that dxSn = 0 in the presence of mixed anomalies for the configurations of the form 
(3.1). It would be interesting to consider more involved geometries for which dxSn might pick 
up the anomaly. Or, perhaps dxSn / 0 if instead of considering the entanglement between 
spatial regions we consider the entanglement between charged degrees of freedom. 
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A Conventions for Wick rotating 

In the majority of this work we carry out computations in Euclidean signature. More precisely, 
the time coordinate in Lorentzian signature is analytically continued to the complex plane 
and computations are carried out along the imaginary time axis. In this appendix we denote 
quantities along the real time axis with a subscript ‘L’ (Lorentzian) and quantities along the 
imaginary time axis with a subscript ‘E’ (Euclidean). In most cases, we work exclusively in 
Euclidean signature and the aforementioned subscripts are omitted. In cases where there is 
a possible ambiguity we specify explicitly which signature metric we are working with. 


- 17 - 


Consider a theory defined on a manifold with Lorentzian signature and time coordinate 
ti- We may analytically continue this time coordinate to a complex one r, 

T = tL-itE, (A.l) 


where ti, and tE are real. We define the Euclidean theory as a restriction of the complexified 
time theory to imaginary time. We denote such a restriction by 


ti —)■ —He ■ (A.2) 

Given an action in Lorentzian signature 



Sl = J dti ^x ^/-gL Cl , 

(A.3) 

we analytically continue it 

so that 

Sl —)• iSe , 

(A.4) 

where 

Se = — J dtE d^~^x Ce , 

(A.5) 

and 

CE{tE) = —Cl{—He) , 

(A.6) 


arranged so that a canonical kinetic term in Ce has opposite sign of that of Cl. 

Given a Lorentzian signature metric gLfiu = V/iu and a boost with rapidity kl, we can 
extend the boost to the complexified metric via k —>■ kl — ike- The angle ke (denoted by 
9 in the main text) will act as a rotation along the complex time coordinate and the sign is 
chosen so that ke corresponds to a counterclockwise rotation in appropriate coordinates. For 
instance, we find that under boosts, zl = x — t is transformed to zl = ZLe.~'^^ in Lorentzian 
signature while its analytic continuation ze = x + He is transformed into ze = To 

emphasize this point let us consider the generator of boosts in Lorentzian signature 

Xidf, = tdx + xdt . (A.7) 

Analytically continuing to t = ti — He and setting = 0 we find 


Xe^im = -i {tEdx - xOt^) ■ (A.8) 

Thus, 

{tEdx - xdtE) ■ (A.9) 

A Ghern-Simons action analytically continued to the imaginary time axis will become 
imaginary. Let us define 




g/ii 

7 ^’ 


^E 




(A.IO) 
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Given (A.5) we find that 



‘‘^E 

(A.ll) 

Thus, given 

Slcs = J df^XL y/-gL ... , 

(A.12) 

we have 

SlCS —^ iSECS , 

(A.13) 

with 

iSecs = j <fxE{-iy/^) {ie^E'"^'^) ■ ■ ■ ’ 

(A.14) 

resulting in 

SECB = -i f d‘^XEy/^e>^j,^-^‘^ ... . 

(A.15) 

The generating function of connected correlators is given by 



ZL[gLf..] = [ DcPlc^^^. 

(A.I6) 


We can analytically continue this generating function to a Euclidean signature metric by 
extending the metric to a complex one and restricting the path integral to contributions from 
its Euclidean component: 


J = e‘^ZE[gE^J.u] ■ 


(A.17) 


The relative factor denotes a possible phase gained by the Euclidean partition function due 
to the change of variables from (pi to (pE- (For instance, the complexified time component of 
the gauge field is imaginary along the imaginary time axis). Such a constant phase will not 
affect any of our computations and we will set it to zero in what follows. Thus, we hnd, for 
example. 


and hence 



. 2 5\nZ 2 6 Wl 

'I _^ _ 

V-9L Sql fliz V~9 l Sql fliz 

2 <51nZ _ 2 6We 

^9E flU y/^ ^9E 


Wl iWE . 


(A.18) 

(A.19) 

(A.20) 


B Anomalous Ward identities and the Chern-Simons term 

In this appendix we work out the Ward identities associated with the anomalous non-conservation 
law for the stress tensor in the presence of anomalies. We present a formal derivation of these 
Ward identities in section B.l using the anomaly inflow mechanism. In section B.2 we rederive 
the Ward identity for gravitational anomalies in three dimensions in an explicit manner and 
in section B.3 we connect this result to a holographic computation of entanglement entropy. 
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B.l Ward identities and anomaly inflow 


We have argued that in the presence of anomalies the coordinate variation of the generating 


function 5^W satisfies 


where 


5^W = i 


IM 




*r^ = 


dies 

5rw 


(B.la) 


(B.lb) 


We will now follow [10] and derive this result explicitly. As opposed to the majority of this 
work, the analysis will be carried out in Lorentzian signature. We then follow the conventions 
of appendix A to analytically continue our results to Euclidean signature. 

The anomaly inflow mechanism, discussed in [50], posits that while / 0 for a theory 
with an anomaly, the variation of the covariant generating function. 


W' = W+ f Ics, (B.2) 

Jd-^M 

satisfies 6 ^W' = 0. Here d~^A4. is a manifold whose boundary is Ai which is where W is 
defined and Iqs is a Chern-Simons form. Put differently, the non-coordinate invariance of W 
is equal to the non-coordinate invariance of the Chern-Simons term, 

S^W = [ Ics . (B.3) 

Ja-^M 

If is a tensor, then the change in under small coordinate transformations is 
given by its Lie derivative 

d^A"^ = r^A"^ . (B.4) 

The Christoffel connection is not a tensor and satsifies 


where we have defined 

and 


+ dVj , 

vj = daX^ , 


Thus, 




= dVj>^ + C^Ics, 


(B.5) 

(B.6) 

(B.7) 


(B.8) 
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where we have written 7cs as a function of the Christoffel connection one-form, the Riemann 
tensor two-form = R°^dx^ and the Abelian gauge held one-form A and its held 
strength F = dA. We use roman indices for components in 

In order to evaluate the Lie derivative of the Chern-Simons form we use 

Ics = ics^ , (B.9) 

where Li is the volume form on Using we have 

J ^x^CS = y = J diy^Ics , (B.IO) 

where denotes the interior product. Thus, under the integral we have 

<5*Ios = -Va^d (1^) + d + i^Ics) . (B.ll) 

Since the Chern-Simons term is invariant under coordinate transformations up to boundary 
terms, we conclude that is closed and we are left with 


d^Ias = d + i^lcsj . (B.12) 

To obtain (B.l) we restrict ourselves to coordinate transformations for which 

X^ = 0, (B.13a) 

5xx“ = 0, (B.13b) 

where T denotes the coordinate in the bulk of d~^A4. Having x"*" = 0 implies that 

/ Cycles = 0, (B.14) 

Ja-^M 

and that I 4 -*- = 0, while (B.13b) implies that Vj_“ = 0. Inserting these relations into (B.12) 
and using (B.3) leads us to 

dxW = -[ (B.15) 

JM 

with defined by (B.lb). 

So far, we have carried out our analysis in Lorentzian signature. To go to Euclidean 
signature we note that each of the terms in 

W'l = Wl +J Ilcs, (B.16) 

is continued to Euclidean signature via 

W'l^iWe, Wl^iWe, I Ilcs^ I Iecs- (B.17) 
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(Where the relation for the Chern-Simons form follows from the second equality in (A.12).) 
Thus, the Wick rotation 

Wl + J Ilcs ^ + j Iecs, (B.18) 

implies 

6 ^WE = iSx J Iecs, (B.19) 

and leads to (B.l) closing our argument. 

The construction which we have used in obtaining (B.l) is a purely theoretical one. The 
anomalous quantum field theory on A4 need not have an extension into d~^A4. Needless to 
say, pion decay into two photons is not considered to be an indication for the existence of extra 
dimensions. For this reason, we have chosen to restrict x according to (B.13) —the theory on 
A4 should not be aware of the precise extension of x from A4 into d~^A4. Indeed, it does 
not seem unlikely that one may add to W boundary counterterms involving the extrinsic 
curvature of Ai from which one may derive (B.l) without the use of (B.13). 

B.2 Gravitational anomalies in two dimensions 

In order to make contact with the works of [23-25, 51, 52] it is useful to consider an explicit 
example. Let us place our quantum field theory on an n-fold cover AAn with d~^AAn being 
an asymptotically AdS space. While our intention is to relate the result of this section with 
computations of entanglement entropy in the context of the AdS/CFT correspondence, one 
can also consider the asymptotically AdS space d~^M.n as a particularly useful extension of 
AAn into an extra dimension with no reference to holography. 

Following the notation of [23] let us denote the Euclidean metric on 9“^A4i_|_e by 

ds^ = + {dx^)‘^) + {g±± + K^x^){dx^)‘^ + e’^'^Ufj,dx^dx^ + ... , (B.20) 

where we have expanded the metric around t = 0, the location of the conical deficit. 

Here // = 0,1. The function (j) serves as a regulator whose derivatives have compact support 
near the tip of the cone and satisfies:^ 

= 4TT6{t, x ^), (B.21) 

which is consistent with R = —dvre with n = 1 + e to linear order in e. 

Following [23] and the convention in (A. 15) let us evaluate the Chern-Simons term 

I Ics = Cgj ^ics = Cg j (dbT^dc + , (B.22) 

on the line element (B.20). In what follows we will often us 6 °“^^ = jy/g with the 

Levi-Civita symbol satisfying e** ^ = 1. As a warmup exercise let us evaluate the term 

®The function 0 satisfying (B.21) is given by <)) = 2 log p where p^ = t^ + ■ A. choice of the regularized 

version of the function is 0 = 2 log \Jp^ o? with a small parameter a [53]. 
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f Iqs on the line element (B.20). After taking into account that (j) has compact support near 
the origin, one finds that the only contribution to terms in f I(js which are linear in e are 
given by 


de y^CsL=o = y 

^ d^xdx^ 

Uc^eP’^dgU, - dg4>5P'^ d,{e^^do,Up) 


= -hj 

^ d^xdx^ Uc^eP'^dgUy, 

(B.23) 

= -2TrCg 

f dx-^eP^dgll,,. 



We would like to see how the result (B.23) changes under a (small) x-^ dependent rotation 
in the t-x^ plane parameterized by 


X“ = 




(B.24) 


That is, we would like to compute -^csIg^Q- Note that (B.24) satishes (B.13a) but not 

(B.13b). We will come back to this point later. The coordinate transformation associated 
with will shift the line element from (B.20) to 


dsg = ds^ + 5gabdx°‘dx^ , 


with 


^Sab ^xSab ■ 

The rotational invariance in the t-x^ plane allows us to write 

5gabdx°‘dx'^ = {5g±±_ + 5K^x^){dx^)‘^ + 5U^^dx^ dx^ + • • 

where, for example, 

5Uy = -69e'^'^U^ + 269'ef^'^a^ + 6^^'^ 69d^U^e'^Pap + ■■■ 


(B.25) 

(B.26) 

(B.27) 

(B.28) 


Since the rotated metric (B.27) has the same structure as the unrotated one (B.20), we 
can use (B.23) directly to obtain 


/ 


Sxde I = -\cg j d^xdx^ [□,/. d,{e'^^d^ 5 U 0 ) 


A short computation yields 


= -Ad9'{x^) + 


(B.29) 


(B.30) 


where • • • denotes terms which vanish when localized at the origin. Therefore, the first term 
on the right hand side of (B.29) will contribute to f Ics, but the second term will not. 
Thus, 





= — -Cg J d^xdx'^ 

=47rCg 69{x'^ = 0), 


(B.31) 
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where in the last equality we have assumed that the entangling region is a semi-infinite line, 
i.e., 66 {x-^ = 0) is single valued. In the case of a finite interval we obtain 

6 ^d, I Icsl^^ = 8TrCg6e. (B.32) 

The results (B.31) and (B.32) may be modified since (B.24) does not respect (B.13b). 
As we discussed, a non-vanishing contribution to f Ics as a result of d±x^ / 0 is a 
remnant of the construction we have been using and may be expected to cancel via boundary 
counterterms associated with the extrinsic curvature of A4„, perhaps similar to the terms 
computed in [24, 51, 54]. Luckily, we do not have to worry about these terms: an analysis 
similar to the one carried out in obtaining (B.23) suggests that 

i>./KL“^ = 0, (B.33) 

implying that no new terms are generated from the violation of (B.13b). 

Using (B.19) we find that the variation of the Euclidean generating function is given by 

6 ^d,Wi+, I = i6^d, I Ics I , (B.34) 

in precise agreement with (2.26) in the limit where n = 1 -|- e. For a finite interval we would 
have obtained 

6 ^d,Wi+,l^^ = 87ricgSd. (B.35) 

in agreement with our field theory expectations. 


B.3 Holography 

We may also use (B.32) to compute the change in entanglement entropy of an interval due 
to a boost as predicted by holography [27, 30, 53]. Consider a bulk action of the form 


'S'bulk — 5'invariant T SqS T 'S'boundary i 


(B.36) 


where ^invariant denotes the fully gauge and diffeomorphism invariant part of the action, 
Scs = -if Ics is the contribution of the Chern-Simons term to the bulk action and ^boundary 
denotes possible boundary terms which may contribute to the on-shell action. Following the 
notation of Dong [53], we decompose Sbuik into a contribution coming from the tip of the 
regularized conical singularity and a contribution coming from outside the tip 


'S'bulk — ‘S'inside T 'S'outside • 


(B.37) 


The prescription of [53] for computing the entanglement entropy can be written in the form 


5a 


^e^inside • 


(B.38) 
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Using (B.32) we find that (for a finite interval and ignoring possible boundary) 


deSA = -dgdeSinside , 

= -ded, (^-i I Ics) , (B.39) 

= Svri Cg , 


in agreement with (2.34). 
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